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ABSTRACT 

We propose that a large fraction of the QSO Lyman hmit absorption systems (LLS) observed at 
high redshift (2 >, 3) originate from gas trapped in small objects, such as minihalos, that form prior to 
reionization. In the absence of a strong UV flux, the gas is predominantly neutral and may form clouds 
with HI column density iVni ^ 10^* cm~^. Due to their high densities and high HI column densities, 
these clouds are not destroyed by the onset of the UV background at a later time. Thus, if not disrupted 
by other processes, such as mergers into larger systems or 'blow away' by supernovae, they will produce 
LLS. We show that the observed number density of LLS at high redshifts can well be reproduced by 
the survived 'minihalos' in hierarchical clustering models such as the standard cold dark matter model. 
The number density of LLS in such a population increases with z even beyond the redshifts accessible 
to current observations and dies off quickly at z ^ 2. This population is distinct from other populations 
because the absorbing systems have small velocity widths and a close to primordial chemical composition. 
The existence of such a population requires that the reionization of the universe occurs late, at z < 20. 

Subject headings: quasars: absorption line - galaxies: formation - cosmology: theory 



1. INTRODUCTION 

Lyman limit absorption systems (LLS), observed as 
breaks in quasar spectra at the rest-frame wavelength 
of the Lyman continuum edge, are produced by gaseous 
clouds optically thick to Lyman continuum photons. 
These systems can be observed over a wide range of red- 
shift, from z ~ to approximately that of the highest- 
redshift quasars. Thus, such systems can be used to probe 
the gas in the Universe at various redshifts. Due to their 
relatively high HI column density (Vhi xL 10"'^^ cm~^), 
these systems are usually thought to arise from photoion- 
ized clouds in gaseous galactic halos (see e.g. Steidel 1990, 
Mo & Miralda-Escude 1996). There is now direct observa- 
tional evidence that some of the LLS are indeed associated 
with galaxies (see e.g. Steidel 1995), at least at redshift 
z < 3. This is not surprising, because any HI disk like that 
of our own Galaxy will produce a LLS when intersecting 
a sightline to a quasar. However, the observed rate of in- 
cidence (i.e. the number of systems per sightline per unit 
redshift) of LLS is large and increases rapidly with red- 
shift (Stengler-Larrea et al. 1995); there is no indication 
of a reversal of this trend even for the LLS with the high- 
est redshift z ~ 4.5 observed so far (Storrie-Lombardi et 
al. 1994). The observed rate of incidence of LLS at z ^ 4 
is about 30 times larger than that expected from the local 
number density of galaxies with absorption cross sections 
given by the optical Holmberg radius. Thus, if LLS at 
z >, 4 were also associated with galaxies, it would mean 
that (1) most galaxies have already formed at that red- 
shift, and (2) the absorption cross sections of galaxies are 
much larger at high redshifts. However, in a hierarchical 
model of structure formation, such as the standard cold 
dark matter (CDM) model and its variants, the number 
density and sizes of galaxies are small at high redshifts. 



because most gas has not yet been assembled into large 
galaxies (Mo & Miralda-Escude 1996). It is therefore un- 
likely that LLS observed at high redshifts are due to galac- 
tic populations. 

It is possible that a large fraction of the observed LLS 
at high redshifts are produced by photoionized clouds in 
small halos. Recent hydrodynamic simulations show that 
low column density absorption systems can be produced 
by gas in filamentary structures that arise from the grav- 
itational collapse of initial density fluctuations (e.g. Cen 
et al. 1994; Zhang, Anninos, & Norman 1995; Hernquist 
et al. 1996; Haehnclt, Steinmetz & Ranch 1996; Miralda- 
Escude et al. 1996). The number densities of LLS given 
by these simulations were, however, found to be lower by 
a factor of about ten than what is observed (e.g. Katz 
et al. 1996). The discrepancy at z ;< 2 may be explained 
by the population of pressure-confined clouds in gaseous 
halos that are not resolved by current simulations (Mo 
& Miralda-Escude 1996). The discrepancy at z >, 3, how- 
ever, was not explained by previous investigations based on 
photoionization models in hierarchical cosmogonies (e.g. 
Gardner et al. 1996). 

It is also possible that the flux of the UV background 
may drop substantially at z > 3 and the HI column den- 
sity of a cloud could be underestimated in models as- 
suming a constant UV flux. However, recent observa- 
tional results based on the proximity effect show that 
the UV flux at the Lyman continuum edge is about 
(1 — 5) X 10~^^ergs~^ cm~^Hz~^ sr~^ at z ~ 4.5 (Williger 
et al. 1994; Giallongo et al. 1996; Lu et al. 1996). This 
value of the UV flux is already comparable to, or even 
higher than, the ones used in most of the model consider- 
ations mentioned above. 

In this paper we suggest the possibility that the discrep- 
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ancy between model predictions and observations of the 
LLS at high redshifts (z >, 3) is a result of the neglect of 
the contribution of niinihalos that collapsed prior to reion- 
ization. The gas in such small systems is predominantly 
neutral and may form clouds with H I column density 
■^Hi ^ 10^* cm~^. Because of their high density and high 
HI column density, these clouds will not be destroyed, or 
ionized by the onset of the UV background at a later time. 
They will therefore produce LLS at lower redshifts, if not 
disrupted by other processes. 

We present our model for the LLS observed at high red- 
shifts and compute their number per unit redshift in Sec- 
tion 2. A brief discussion of our model is given in Section 3. 

2. A NEW POPULATION OF LYMAN LIMIT SYSTEMS 

For simplicity we assume the Universe to be flat with the 
cosmic density parameter ilo = 1. The Hubble constant is 
written as Ha = 100/ikms^^ Mpc^^. 

2.1. Physical properties 

Prior to the epoch of reionization the physical proper- 
ties of cosmological structures can be estimated from the 
spherical collapse model (see e.g. Padmanabhan, 1993). 
The gas in a collapsing dark matter perturbation is shock 
heated to the virial temperature 
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where Vc is the three dimensional virial velocity of the halo, 
mp the proton mass, and fi the mean mass per particle. 

The total H number density in a dark matter halo that 
collapses at a redshift, z, with an overdensity 6 at virial- 
ization reads 
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where Pcrit denotes the closure density. Assuming that 
dark halos are singular isothermal spheres with density 
profile p{r) ex r~^, the virial radius r„, defined as the ra- 
dius within which the mean mass overdensity is IStt^, can 
be written: 
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where t = {2/3Ho){l + z) ■^/^ is the cosmic time at z. 

For halos with T„„. < 10'' K (Vc<15kms-i) the 
gaseous component cannot be coUisonally ionized and will 
remain neutral. Due to the high collapse redshifts (z > 5) 
considered here, the recombination time scale in miniha- 
los is shorter than the Hubble time, and the gas in clouds 
with Vc ^ 15kms~^ before reionization will also be neu- 
tral shortly after their formation. As a result, the column 
density of neutral hydrogen atoms, Nm, is given by 
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The corresponding total mass of such a halo is 
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It is clear that the optical depth at the photoionization 
threshold, 
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exceeds unity for a wide range of collapse redshifts and 
virial velocities. 

To determine what happens to the clouds once the 
UV background switches on we analyze the equilibrium 
equation of the ionized fraction for a pure hydrogen gas, 
X = rie/n, where rig and n denote the free electron and 
total baryonic number density, respectively. This fraction 
is given by 



= kphn{l - x) 
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where kph — 47rJ (J//ii/) exp[— iVHio'(i^)](T(i/)rfi^ and 
kpec denote the rate coefficients for photoionization of neu- 
tral hydrogen and radiative recombination to H I, respec- 
tively. Motivated by the indirect measurements from the 
proximity effect (see Giallongo et al. 1996, and references 
therein) we assume the background UV flux to be given 
by a powerlaw 
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where h is Planck's constant and hvTH = IRyd = 13.6 eV 
is the H I ionization threshold. In the optical thick limit 
the photoionization rate coefficient can be approximated 
by a power law kph ~ J2ikp/^[NH{1 — x)]~^, where Nh 
is the total column density of hydrogen nuclei. The ex- 
ponent P has a non-trivial dependence on a and is deter- 
mined numerically to vary from 1.35 to 1.69 for 1 < a < 2. 
Inserting this approximate expression for the photoioniza- 
tion rate into the equilibrium equation for the free electron 
fraction [equation (7)] yields the implicit solution for x 
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which for (3 — reduces to the more familiar optical 
thin case. We adopt a = 1.8 for the spectral index, 
with which the photoionization rate can be approximated 
by kph « 4.1 X 10-^^J2i{Nm/lO^^ cm-^)-/^ g-^ with 
/3= 1.6. 

The right hand side of equation (9) has a maximum at 
X = Xmax ~ 0.8 with a value of f{xmax) ~ 0.25. The exis- 
tence of this maximum shows that ionization equilibrium 

cannot be achieved for J2ikpf^N(^^ /{nkpec) > J{xmax)- 
In this case photoionization overcomes recombination and 
ionizes the gas until a new, highly ionized, equilibrium 



state is reached at column densities ^ 10^^ cm^^. As 
discussed in Miicket and Kates (1997), the evolution of 
such systems from one equilibrium state to another can 
be rapid. The recombination rate'^ suitable for an equilib- 
rium temperature of 10'* K is fc/jec ~ 2.6 x 10~^^ cm"^ s~* 
(Ferland et al. 1992). Thus, for a given J21, the maximum 
of f{x) can be used to derive a critical total hydrogen 
column density NfJ"^'^. Note that this column density is 
calculated at the edge of the cloud and is approximately 
equal to ryirUvir [see equation (2) and equation (3)]. We 
can therefore define a minimum collapse redshift, Zcrit, for 
halos of a given virial velocity, so that they will remain 
optically thick if they form at z > Zcrit- This minimum 
collapse redshift can be written as 
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As we will show below, minihalos (with Vc ^ lOkms^^) in 
current hierarchical models have typical collapse redshifts 
higher than this critical value; they can therefore retain 
their high initial H I column densities. 

To check the robustness of the above derivation for 
the critical collapse redshift wc have integrated the time- 
dependent chemistry and cooling model presented in Abel 
et al. (1997a). We found good agreement between the nu- 
merical results and the analytic expression of the critical 
redshift given by equation (10). 

The density of gas contained in minihalos considered 
here is quite high due to their high collapse redshifts. 
The recombination time scale for the gas is short, about 
0.1 Myr, which makes them stable against photoioniza- 
tion from internal stellar sources. On the other hand, 
only a few supernova explosions may be able to blow out 
most of the gas from the halos (see Ciardi & Ferrara, 
1997, and references therein) if star formation can hap- 
pen in them. However, systems that form in halos with 
Vc <, 15kms^^ at z ^ 15 are not be able to cool since H2 
formation is inefficient in these halos at such redshifts (see 
Tegmark et al. 1997, and Abel et al. 1997b). Also the 
UV flux of the first structures in the universe may lead to 
H2 photodissociation prior to H I reionization (Haiman, 
Rees and Loeb 1997), inhibiting line cooling by H2 . Fur- 
thermore H I line cooling is not efficient due to the low 
virial temperature {^ 10* K). Thus these systems may 
not form stars and hence are not subject to "blow away" . 
The gas can only be gravitationally confined in the dark 
matter halo if the typical particle velocity is less than 
the escape velocity of the halo. The mean particle ve- 
locity of hydrogen atoms in a Maxwellian distribution is 
Vh{T) = {SkT/nipy^'^ for a neutral gas, and the escape 
velocity is given roughly by Vesc ~ V^Vc- Hence only 
systems with Vc >, Vh{10'^ K)/-\/3 « 9kms~* are stable 
against evaporation. These considerations suggest that 
minihalos with Vc — 10-20 km s~* are most likely to pro- 
duce LLS. 




Fig. 1. — The numfeer orLLS per 'unit* redshift. The curves show 
the prediction of the Standard CDM model with erg = 0.7. Results 
are given for minihalos in three Vc bins. The proposed minihalo 
LLSs provide a significant fraction of all LLS only at redshifts 2^2. 
The observational data is taken from Stengler-Larrea et al. (1995) 
(filled circles) and Storrie-Lombardi et al. (1994) (filled squares). 

2.2. The predicted number density of LLS 

The typical collapse redshift of halos with 
Vc ^ SOkms^* in current cosmogonies such as the stan- 
dard CDM model exceeds the critical redshift given in 
equation (10). Such small halos are also abundant at high 
redshift in these models. It is therefore possible that mini- 
halos contribute a substantial part of the total number of 
LLS at high redshifts. Here we examine this possibility. 

We assume that minihalos which form at z > Zcrit and 
do not merge into larger systems by a redshift Zs will con- 
tribute to the absorption cross section at Zs- The con- 
ditional probability, p{M2, Z2\Mi, zi)dM2, that a halo of 
mass Ml selected at redshift zi will have merged to form a 
halo of mass between M2 and M2 -I- dM2 at a later redshift 
Z2, can be calculated from the extended Press-Schechter 
formalism (Bond et al. 1991; Bower 1991). As in Lacey & 
Cole (1994), halos are assumed to be destroyed by redshift 
Zs , if they have merged into halos of twice their initial mass 
by that redshift. The fraction of minihalos that form at 
zi and survive merging until Zg, fs, can then be obtained 
from 
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The absorption cross section of a minihalo is assumed to 
be spherical with radius given by equation (3). The re- 
sult for a standard CDM model normalized to as = 0.7 
is given for three different bins of virial velocities in Fig- 
ure 1. In the calculation zi is chosen as the redshift at 
which the comoving number density of halos prediced by 
the Press-Schechter formalism peaks. This corresponds 
to the collapse redshift of a ~ 1.5(7 peak. In all cases, 

Zl ^ Zcrit • 

Clearly, the total number of systems with 
10kms~* <, Vc ;S 15kms~* is sufficient to explain the 
observed number of LLSs at z >, 3. This is consistent with 
the observational evidence that LLS at low redshift are as- 
sociated with galaxies (see Steidel et al. 1996). As shown 
in Mo & Miralda-Escude (1996), clouds pressure-confined 
in galactic halos can indeed give a sufficiently large cross 
section to explain the number density of LLS at z <, 2. 
Note that for current low-f2 cosmogonic models having 
smaller power on small scales than the standard CDM 
model, the number density of minihalos may be reduced. 
However, the predicted number of absorption systems per 
unit redshift may not be reduced substantially, because 
it depends also on dl/dz (the proper distance per unit 
redshift) and the virial radius, both are larger in a low- 
rt model. The dependence on cosmological parameters 



^Note that the Case B recombination rate is used for the optical thick case considered here. 



cannot be examined in detail until more detailed hydro- 
dynamical calculations are able to determine the physical 
sizes of the absorbers. 

It is important to realize the difference between our 
model and the minihalo model proposed by Rees (1986) 
for Lyman alpha forest systems. In the model of Rees, 
minihalos are optically thin, highly ionized and have an HI 
column density of ^ lO'^^ cm^^, characteristic for Lyman 
alpha forest systems. This assumption is correct for mini- 
halos that collapse after reionization from an already ion- 
ized intergalactic medium at redshifts z < Zcrit [see equa- 
tion (10)]. 

3. DISCUSSION 

As shown in Thoul and Weinberg (1996), after reion- 
ization the collapse of gas can only happen in halos with 
K;;^30kms^^, and so minihalos (with Vc/^SOkms^^) 
that form after reionization are not expected to give a 
significant contribution to the LLS. In contrast, prior to 
reionization the collapse of gas can even happen in halos 
with Vc as small as about 0.1 kms^^(l -I- Zcoii)- Our model 
for the LLS envisages a hierarchical cosmogony where a 
large number of minihalos (with Vc ^ 3Gkms~^) can form 
prior to the epoch of reionization. We show that the gas 
trapped in such minihalos is, due to its high density, able 
to withstand photoionization by the UV background set 
on at a later time and reaches thermal and chemical equi- 
librium at a temperature ~ 10^ K. By including these 
systems, the discrepancy in the number of LLS between 
the observational result and earlier theoretical predictions 
based on CDM models can be reconciled. The model pre- 
dicts the number density of LLS to keep increasing towards 
higher redshifts. 

The minihalo population is assumed to form at a time 
when the gas in the IGM is largely primordial. As a re- 
sult, the absorbing gas in these systems should also have a 
close to primordial composition. Although later accretion 



of gas from the enriched IGM may increase the metallicity, 
we expect the metallicity of such LLS to be low. 

Minihalos with lOkms^^ <, Vc ^ 15kms^^ are ex- 
pected to dominate the total cross section of LLS at high 
redshifts, because (i) they can form in the low-pressure 
IGM prior to reionization, (2) they are stable against evap- 
oration, (3) they cannot cool by II2 or H I to form stars 
and so are not subject to supernova explosion or inter- 
nal photoionization, and (4) their total cross section is 
large in hierarchical models. Thus, the minihalo popu- 
lation of the LLS should be observed to have a narrow 
range of velocity widths (with H I Doppler parameters 
of <,20kms~^). This can be tested by high-resolution 
spectroscopy of LLS at high redshifts. The two LLS at 
redshifts 3.32 and 2.80 studied by Tytler et al. (1996) us- 
ing the HIRES spectrograph on the KECK telescope have 
Doppler parameters and velocity offsets between individ- 
ual components less than 2Gkms~^. This is consistent 
with our model. Obviously more high-resolution data at 
high redshifts are needed to constrain the fraction of LLS 
in the minihalo population. The sizes for the minihalo ab- 
sorbers are typically about 1 kpc, which may be tested by 
the spectra of gravitationally-lensed quasar pairs at high 
redshifts. 

Much theoretical work remains to be done to model the 
hydrodynamic and radiative-transfer processes to obtain 
more accurate predictions for the properties of the mini- 
halo Lyman limit systems proposed here. 
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